Cold Fermi atomic gases in a pumped optical resonator 
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We study systems of fully polarized ultracold atomic gases obeying Fermi statistics. The atomic 
transition interacts dispersively with a mode of a standing-wave cavity, which is coherently pumped 
by a laser. In this setup, the intensity of the intracavity field is determined by the refractive index of 
the atomic medium, and thus by the atomic density distribution. Vice versa, the density distribution 
of the atom is determined by the cavity field potential, whose depth is proportional to the intracavity 
field amplitude. In this work we show that this nonlinearity leads to an instability in the intracavity 
intensity that differs substantially from dispersive optical bistability, as this effect is already present 
in the regime, where the atomic dipole is proportional to the cavity field. Such instability is driven 
by the matter waves fluctuations and exhibits a peculiar dependence on the fluctuations in the 
atomic density distribution. 

PACS numbers; 42.50.Pq,42.50.Nn,42.50.Wk,71.10.Fd 
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I. INTRODUCTION 

Cavity Quantum Electrodynamics (CQED) is one of 
the most rapidly developing areas of modern quantum 
optics and atomic molecular, and optical physics P, Q. 
In its early days optical instabilities, with the paradigm 
example of optical bistability were in the center of inter- 
est [3, Ijl . Optical bistability emerges from the nonlinear 
response of the dipolar transitions of atoms to the cavity 
field, to which they couple. Theoretically, such dynam- 
ics are well described by mean field theories formulated 
for macroscopic variables, such as the cavity field inten- 
sity and the total atomic polarization [1, [g, 01 ■ Fluctu- 
ations of the polarization, in particular in the quantum 
regime, may play an important role in these processes, 
as pointed out for instance in the seminal studies of the 
late Dan Walls and collaborators (see Q and references 
therein). Typically, however, these dynamics are charac- 
terized by few macroscopic quantum degrees of freedom, 
where quantum fluctuations are small, except at the in- 
stability points. Nonlinearity at the microscopic level 
has since then been reached by high- finesse resonators, 
where the strong coupling regime between a single atom 
and a single photon have been realized in several mile- 
stones experiments [lo| . In this regime the quantum 
dynamics of few microscopic quantum degrees of freedom 
dominates the behavior of the system, and quantum fluc- 
tuations are typically large. 

Recently experimental progress has allowed one to 
strongly couple cold atomic gases with the electromag- 
netic field mode of a resonator. This line of research, 
which we would like to term cavity QED of many-body 
systems, was stimulated by the observation on the inter- 
dependence between intracavity field and atomic motion 
via of the mechanical effects of light ,11, 12]. This was 
confirmed by seminal experiments, which demonstrated 
cavity cooling of atoms [13l| and selforganization of atoms 



in regular patterns inside pumped cavities [Tj, [T^ . Most 
recently, the coherent coupling between a Bose-Einstein 
condensate and a cavity field has been experimentally 
demonstrated [l^, [l3, ES]. Moreover, the nonlinear in- 
teraction between ultracold atoms and the field of a cav- 
ity mode had been experimentally observed [l9j . In this 
new ultracold regime, theoretical work on bosonic atomic 
gases showed that the matter- wave quantum fluctuations 
and the spatial mode variance give rise to additional non- 
linear atom- field effects [1^, leading to features of the 
phase diagram which suggest novel bistability phenom- 
ena [Ullla,!!!]. 

In this work, we consider a gas of ultracold and polar- 
ized fermionic atoms, whose dipolar transition interacts 
dispersively with the mode of a standing-wave resonator. 
In the regime in which the atomic polarization is linear 
in the field amplitude, and hence where there is no usual 
optical bistability, we encounter instabilities of the intra- 
cavity field which originate from the coupling with the 
atomic motion. The theoretical model we apply is based 
on a second-quantized approach for the description of the 
atomic and of the field degrees of freedom, in the regime 
in which the atomic internal degrees of freedom are elim- 
inated from the equations in second order perturbation 
theory. These dynamics are studied by numerically inte- 
grating the Heisenberg equations of motion. The results 
show that the instabilities of the intracavity field inten- 
sity are driven by the pump strength and by the atomic 
density, and indirectly depend on the degree of localiza- 
tion of the atoms at the minima of the cavity potential. 

This article is organized as follows. In Section |TT] we 
theoretically derive the model, in Sec. IIIII we solve nu- 
merically the coupled Heisenberg equations of motion for 
field and atoms for a fixed number of atoms and thor- 
oughly discuss the validity of the considered approxima- 
tions and experimental parameters. The conclusions and 
outlook are presented in Sec. IIVI 
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II. THEORETICAL DERIVATION OF THE 
MODEL 

In this section we derive the theoretical model, which 
is at the basis of the studies presented in Sec. IIIII In 
order to introduce the many-body theory, which is the 
starting point of our study, let us first review the Hamil- 
tonian of a single two-level atom of mass m interacting 
with a cavity field which is pumped coherently by an ex- 
ternal, monochromatic field. We denote by tr = |5)(e| 
and (7^ the atomic lowering and raising operators of the 
dipole transition at frequency Ua with ground and ex- 
cited states \g) , |e), and by a and the annihilation and 
creation operators of a cavity photon at energy hujc- We 
further assume that the atomic motion is along the cav- 
ity axis X, while the radial degrees of freedom are frozen 
out, and denote by the operators P and X the atomic 
centre-of-mass momentum and position. In the dipole 
and rotating wave approximation, the Schrodinger equa- 
tion of atomic and cavity degrees of freedom is governed 
by the Hamiltonian [2l[ 



H 



JC 



2m 



—ififj {a — a^) . 

(1) 

which is here given in the reference frame rotating at 
the pump frequency LOp. The parameters Aq = ujp — uja 
and Ac = ujp ~ lOc are the detunings of the pump from 
atomic and cavity mode frequencies, respectively, and 
f) is the pump amplitude. The dipole-cavity mode cou- 
pling has position-dependent strength g{X) — go cos{qX) 
with q the cavity-mode wave number. For large detun- 
ing |Aa| the internal atomic level can be eliminated from 
the equations of cavity and atomic external degrees of 
freedom The effective Hamiltonian reads 

H = \- hUna'^acos^iqX) - hAca''a 

2m (2) 

—ihrj (a — dt) , 

with Uq = gg/Aa, and where we have considered temper- 
atures for which the atomic motion is frozen on the time 
scale in which the internal degrees of freedom appreciably 
evolve. 

For later convenience we introduce the dimcnsionless 
variables 



optical bistability. Dispersive optical bistability, in fact, 
would appear when the next order in the perturbative ex- 
pansion is included, giving rise to a nonlinear (Kerr) term 
of the form d'^d^dd [24, 25]. In this work we deal with the 
parameter ranges in which such Kerr nonlinearity can be 
neglected, and hence classical optical bistability does not 
appear. Moreover, we consider the regime in which the 
mechanical effects of photon- atom interactions are rele- 
vant. In this range, the nonlinearity in the dynamics of 
Hamiltonian ^ is solely due to the coupling between 
the cavity field and the atomic quantum center-of-mass 
variables. 



A. Many-body Hamiltonian 

We now introduce the system which we investigate in 
this paper, namely a gas of A'^ fully polarized fermionic 
atoms interacting dispersively with the mode of the res- 
onator. Denoting by ^'(x) the atomic field operators 
obeying Fermi commutation rules, the rescaled Hamil- 
tonian in second quantization reads 

H ~ Acd^d ~iri{a~ a)) (4) 

-I- j dx(^i!\x)-^'^{x) + UoCOS^{x)a)i>\x)ii{x)a 

where we have taken care of the proper ordering between 
atomic and field operators, which is here normal (see 
also [12). 

Let us now consider that the number of photons is fixed 
to the value n. Then, after tracing out the photonic 
degrees of freedom from Hamiltonian (j4|), the resulting 
Hamiltonian describes the motion of N Fermions in the 
optical lattice of the cavity field, whose depth is propor- 
tional to the number of photons n. Hence, the degree of 
localization of the atoms at the minima of the potential 
wells depends on n, and so will the Wannicr functions in 
the tight-binding limit. We now assume that the atoms 
are well localized at the minima of the potential, and 
expand the atomic field operators in the tight-binding 
limit [26j in the corresponding Wannier functions for the 
lowest Bloch band, 



K 



^{x) = ^ fiWn{x - Xi) 



(5) 



p = X = qX, t = tUJr 

nq 



jj Uo Ac V 

Uq = , Ac = , 1] = . 



(3) 



where Wr — is the recoil frequency. 

We remark that when neglecting the coupling between 
cavity and motion, hence in the classical limit, Hamilto- 
nian in Eq. Q does not give rise to instabilities like in 



where fi is the annihilation operator of a fermion at site 
i, Wfiix — Xi) is the Wannier function localized at site i, 
and K is the number of lattice sites. The subscript n 
indicates the operator valued dependence of the Wannicr 



functions on the number of photons n 



Such an 



expansion, considering only the lowest band, is known to 
be justified for large or moderate lattize depths, which we 
will discuss in more detail later on. Using the Wannier 
expansion within the single-band approximation, Eq. 
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we obtain from Eq. (|4]) the Hamiltonian 

K 

(6) 



-AcCi^a — irj {a — a^) 



where 



Eij{n) I dxwlix- Xi) [--^] Wn{x- Xj), 



j{n) ~ J dx w^{x — Xi) cos^{x)wn{x — Xj), 



(7) 



are the coupHng parameters which depend on the number 
of photons of the cavity field. In the tight-binding ap- 
proximation we keep only on-site and nearest-neighbour 
interactions in Eq. ([6]), and denote the relevant coupling 
parameters by £^ = En, J = Ju, Ei = Ei^i+i, and 
Ji = Ji^i+i- This is legitimate in the single band regime. 
We impose periodic boundary conditions and use the rep- 
resentation Cfc = Xlj/i^'*^"'; where (c^) annihilates 
(creates) a fermion at the rescaled quasi-momentum k, 
with k in the interval [—1,1] and step Sk = 2/K. In 
this representation, the number operator N = J^j fjfj 

and the hopping operator B = fjfj+i + H-C- get the 
diagonal forms 



B = '2j2k cos(fc7r)c^Cfc, 



(8) 



Hence, the new Hamiltonian takes the form 

H^-Aca''a-ir]{a-a^) + niih)+H2in)a^a (9) 

with 

Hi{n) = E + cos(A:7r) 



n2{n) ^UoJ2 [Jif^-'^) + Ji(n- l)2cos(A:7r) 

(10) 

Here, we used the relation J-{h)d — d!F{n — 1) for any 
function !F{z) which is analytic in the scalar variable z. 



B. Nonlinearities in the coupling parameters 

Equation Q describes dynamics, which couple in a 
nonlinear way cavity and atomic degrees of freedom. In 
particular, the coupling parameters E, Ei, J, Ji depend 
on the number of photons, since they are integrals of 
Wannier functions, which themselves depend on the in- 
tensity of the cavity field. In order to better understand 



the character of this dependence on the intracavity pho- 
ton number we evaluate their explicit form by replacing 
the Wannier functions in Eqs. ([7|) by Gaussians functions. 



Wn{x - Xi) W w'^{x - Xi) 



1 



(11) 



where a = \V\~^^'^ and V = Uqu is the potential depth 
at a fixed number of photons n. Moreover, we mod- 
ify the Gaussian functions by imposing the condition 
/ dxw'?{x — Xi)w^{x—Xj) = Sij, which allows us to avoid 
non-physical contributions. This ansatz is usually valid 
within the tight-binding approximation as will be fur- 
ther discussed in Sec. IHIDl and has been checked in [2^ . 
Using these modified Gaussian functions we solve the in- 
tegrals in Eq. ^ and obtain 

y 



J = - (l-se- 
2 ^ 



(12) 



Ji 



s—e ■If e 
2 



where we have introduced y = = \V\'^^^'^ and 
s — sign(Aa). Relaxing the orthogonality condition 
one would find that the amplitude | Ji | is different for 
red or blue detuning. The sign dependence s in J and 
Ji arisses from the fact that the Wannier functions are 
centered at different positions in the two cases. From 
these expressions we find |Ji| < 18 J within the valid- 
ity of the tight-binding approximation, which according 
to Ref. [1^ is given by y < 0.5 — 1. Moreover, the 
contribution from ^-neighbours can be estimated to be 
\Ji/Ji \ = exp [ — (^^ — l)f^] ^ 1, which justifies neglect- 
ing next nearest-neighbour terms and so forth in Eqs. ([9]). 



III. DYNAMICS AND STEADY STATE OF THE 
CAVITY FIELD 

In this section we study analytically and numerically 
the dynamics and steady state of the cavity field, fo- 
cussing on instabilities due to the interplay between field 
and matter wave fluctuations at a fixed number of atoms 
N. We consider the Heisenberg equations of motion for 
the field operators, 

a — ~i[d,T-li{'h)] — i[d,H2{n)]n 

—iH2{n)d + iAcd — nd + r] + £_a, (13) 



r] [d + d'^) — 2Kh + 



(14) 



and similarly for a^. These are coupled to the equa- 
tions for the atomic field operators, and where S,i are 
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the quantum noise operators, such that (^i(i)) = and 
{i^it)ij{t')) = 2KS{t - t') for i = a,n. Such fluctuations 
vary on time scales which are faster than 1/k, and will 
be neglected when taking the mean values on the coarse- 
grained time scales, at which we analyze the dynamics in 
what follows. The numerics are carried out by first evalu- 
ating the commutators and then replacing the operators 
by c-numbers, a — > a, — > a* and n ^ n — and 
then solving the differential equations using the regular 
Runge-Kutta method. As initial conditions we thus use 
n{t = 0) = no, a = y^ri^ and a* = y/n^ for some real 
uq. Note that the commutators are carried out exactly, 
i.e. [a,g{n)] = a {gijij — g{n — 1)), and the c-number 
replacement is done afterwards. Truncating the system 
dynamics to only three equations and dropping the op- 
erator structure should be justified for the open system 
we consider. 



A. Intracavity field intensity 

The intracavity intensity is found from the mean value 
of photons. 



n{t) = {h{t)). 



(15) 



This average is taken over field and atomic degrees of 
freedom, assuming that the cavity field variables relax to 
a stationary value on a faster time scale than the time 
scale of the evolution of the atomic variables, hence for 
K ^ kBT/h, with K the cavity decay rate, ks Boltz- 
mann's constant and T the temperature. In this regime, 
a closed expression for the photon number operator can 
be found in the limit in which the number of photons is 
sufficiently large (In fact, the analytic expression gives to 
a good approximation also the steady state solution for 
weak fields.). In this regime we approximate 



(16) 



where is any analytical function and d indicates deriva- 
tive with respect to photon number, and from Eq. (jlOp 
we find 



cick. 



'H2{n) w [/o^ \{J -dJ) + (Ji - 9Ji)2cos(fc7r) 

k 

(17) 

Using similar relations for the commutators in Eq. p4p . 
we find an analytical expression that the mean photon 
number has to fulfill. 



where 



£,{n) =n2 + dn2 + dUi + dn2n. 



(18) 



(19) 



In order to evaluate Eq. (1151) '^^ assume that the atomic 
state is the ground state of the potential, whose depth 



is determined by the time-dependent number of photons 
n(t). The ground state is then evaluated assuming that 
the tight-binding and the single-band approximations are 
valid during the evolution of the field. This treatment is 
hence valid in specific parameter regimes, namely when 
the rate of variation of the photon number, giving the 
height of the cavity potential, remains smaller than the 
energy gap between the lowest and the first band, and 
when the mean number of photons n warrants sufficient 
localization of the atoms, according to the relation 



1 



< 1. 



We have numerically checked that for y < 0.5 (which 
coincides with the tight-binding and Gaussian approx- 
imations), the first band gap is larger than the recoil 
energy. Thus, remaining in the validity regimes for the 
tight-binding and Gaussian approximations (y < 0.5), 
also warrants no non-adiabatic coupling of excited bands. 
Further, the large band gap makes the system resistant 
to nonzero temperature excitations. 



B. Many-body atomic state 

In order to calculate any quantities we need to specify 
the particular state of the atoms. In the regime in which 
we can assume that the cavity field variables relax to a 
stationary value on a faster time scale than the time scale 
of the evolution of the atomic variables, we can eliminate 
the photon variables from the atomic Heisenberg equa- 
tions. In doing this we follow the procedure discussed 
in p^ . now applied to fermionic atoms. Hence finding 
that the atomic dynamics are determined by an effective 
Hamiltonian of the form 



n,s = EN + fiN) 



El 



where 



/(iV) = — arctan 

K 



CiN) 



B, (20) 



CiN) = A, - UoJoN 



(21) 



and the coupling terms given in p2p . This Hamiltonian 
is derived at leading order in an expansion in 1/A^, as- 
suming a scaling such that the number of lattice sites K 
is proportional to the number of atoms N where we as- 
sume a fixed number of atoms N. The coefficients of this 
Hamiltonian depend on the atomic variables through the 
nonlinear coupling with the cavity field [2l|, [S^l . When 
they are independent of the atomic variables, for instance 
in the case of Fermi gases in deep optical lattices in free 
space, the ground state of this Hamiltonian is the usual 
Fermi sea, where the value minimized. Thus, 
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the many body ground state is found by fiUing with one 
atom all states with quasi-momentum k up to the Fermi 
quasi- momentum kp, assuming that kp is smaller than 
the edge of the Brillouin zone. This simple procedure 
cannot be applied when the nonlinearity in the Hamilto- 
nian gives dispersion curves that may exhibit local min- 
ima, which can be at values of k different from k = 
(at the center of the Brillouin zone). Finding the dis- 
tribution that gives the smallest energy is in general a 
linear programming (LP) optimization problem with two 
constrains: (i) at most one particle can occupy each site 
and (ii) the total number of particles equals N. How- 
ever, since the number of particles per site is a discrete 
variable, = 0, 1, standard LP techniques are not appli- 
cable. We use instead a variational approach where we 
assume that the usual Fermi sea is the ground state, and 
check its stability against perturbations. Such check is 
made by randomly extracting a single fermion from the 
Fermi sea and putting it in a random state outside the 
sea. We then compare the energies, calculated using Hcs 
in Eq. (PO)) . for the Fermi sea and the perturbed state. 
Repeating this procedure sufhciently many times, we find 
that the energy of the Fermi sea (according to our ansatz) 
is usually lower than the perturbed one. In some cases, 
however, the perturbed state is energetically favorable. 
On the other hand, the difference in energy between the 
two states is extremely small for the parameter regimes of 
interest and hence does not affect the general structure of 
our results. We have also compared the energy of a fully 
random distribution with the one of the Fermi sea and in 
all examples we considered found the usual Fermi sea to 
have a lower energy. Still, these situations are potentially 
of great interest, since they may indicate novel quantum 
phase transitions due to the change of geometry, or even 
topology of the Fermi surface, similar to what happens 
in graphene (cf. [l^,!!!!)- These questions go beyond the 
focus of the present paper, and will be studied elsewhere. 

C. Quantum optical bistability 

We now evaluate Eq. assuming that the atoms are 
in the state described by the usual Fermi sea, and find 
a nonlinear equation for n, which may exhibit multiple 
solutions. We remind that such nonlinearity is not due to 
the nonlinear coupling between atomic dipoles and fields, 
but instead originates from the quantum fluctuations of 
the atoms at the minima of the confining potential. In 
order to highlight this behaviour, we express n = I/Uqh^ , 
using the relation between the depth of the potential and 
the width of the Gaussian functions. From Eq. (fT5|) . us- 
ing the approximation [a^g{n)\ « ^^^o, = dg{n)a in 
Eq. (fT9|) . we find the nonlinear equation which deter- 
mines y, 

- |a, - ^ (/i(y) + Be-^h{y)) |' 

(22) 



where B = 2^^^-^ cos(fci7r)/iV is the mean value of the 
operator (B) /N, taken on the Fermi sea of atoms, and 

hiy) = l + sy-se-yil + y/2) (23) 

Here we neglected small terms such as the derivative of 
the El, Ji-terms, which indeed has been verified to only 
give minimal corrections. 

Figure [T] displays the set of values (n, N) for which 
Eq. ((22|) is satisfied. The parameter regime is restricted 
to y < 0.5, where our treatment is valid. In particular, in 
Fig.[l](a) the detuning > and the atoms are trapped 
at the nodes of the cavity potential, where for maximum 
localization, y 0, the coupling with the resonator van- 
ishes. Nonlinearity here arises due to the finite size of 
the fiuctuations. In Fig. [Ijb), instead, Aa < 0, and the 
atoms are confined at the antinodes, where the coupling 
to the field is maximum. For large atom numbers the 
effect of small position fluctuations at these points gives 
a small correction to the cavity intensity, and the sys- 
tem shows less sensitivity on the external parameters, 
see also f2ll.l2^. 
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FIG. 1: Set of parameters (n, N) for which Eq. (f22|l is sat- 
isfied. The dimensionless parameters are n — 1, and (a): 
Uo = 10, »7 = 10, Ac = 10 and (b): Uq = -1, rj = 30, 
Ac — —20. In both plots we assumed K = 50 sites of the cav- 
ity potential. The curves are shown for the values of (n, A'^) 
such that y < 0.5, when the tight-binding approximation is 
valid. 

In order to check the stability of the solutions, we nu- 
merically solve the Heisenberg equations in Eqs. (jl3p - 
(I14p for various initial values of n at a fixed number of 
atoms N, as previously explained. Figure [5] displays the 
evolution of some different values of the mean number 
of photons n{t). In particular. Figs. [Ha)-(c) have been 
evaluated in the same parameter regime of Fig. [T](a) and 
for a particular choice of the number of atoms N. In 
Fig. [5] (a) and (b) it is visible that the mean value ex- 
hibits two possible asymptotic solutions, which depend 



6 



on the initial value, while in Fig. [2] (c) there is only one 
asymptotic value of n. Comparing with the results in 
Fig. [Ija), we find that only the solutions of the lower 
branch of the curve {n, N) in that figure are recovered, 
while the upper branch is unstable. The additional solu- 
tions with small asymptotic n in fig. [2] (a) and (b) corre- 
spond to y > 0.5 and are thus beyond the validity regime 
of our approximations and hence not visible in Fig. [1] 
(a). The lower two plots. Figs. [2ld) and (e), should be 
compared with Fig. [TJb). Here, one can observe that a 
small change in atomic number, = 20 or TV = 18, re- 
sults in intracavity field intensities differing by a factor of 
about 100. We again pointed out that by calculating the 
Heisenberg equations of motion we use the exact relation 
[a, g(n)] ~ {g{n) — g(h — 1)) a, and do not impose any 
approximations involving derivatives as in (jl6p . Still, the 
asymptotic solutions of the numerical integration agree 
well with the ones obtained from ([H]), which verifies the 
trustworthiness of the steady state results [H]) . 




t 



FIG. 2: Mean number of photons n as a function of time, 
obtained by numerically integrating the Heisenberg equations 
in Eqs. (|13p - (|14|l for different initial values n(0) and for a 
fixed number of fermionic atoms TV. In (a)-(c) the parameters 
correspond to the ones of Fig. [lla), in (d),(e) to the ones of 
Fig. [Tfb). In the calculations we assumed that the atoms are 
in the ground state of the cavity potential, and that the tight- 
binding regime and the single-band approximation are valid. 
This is fulfilled for n > 0.4 in (a)-(c), and for n > 4 in (d),(e). 
Once the curves reach below these limiting values, they lie 
outside the validity regime of our approximations. 

We finally study the instability of the mean number 
of photons n, and hence the cavity potential depth, as 
the pumping amplitude is varied. Typical examples are 
shown in Fig. [31 where the solutions for n, N of Eq. ([221) 
are plotted as a function of the pump amplitude rj. The 
dashed middle branch represents an unstable solution of 
the equations. Note that the lower branch is almost iden- 
tical to n = 0. We see that n exhibits jumps at criti- 
cal values of the pump 77, corresponding to an instabil- 
ity due to the nonlinear behaviour, to values in which 



our approximations are invalid. Such regime is in the 
grey-shaded zone. Here, our treatment is invalid and we 
cannot make definite predictions. We conjecture that at 
these points the fermions are populating higher Bloch 
bands or forming a Fermi liquid. The study of this tran- 
sition will be explored in detail in future works. 
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FIG. 3: Mean number of photons n, evaluated in Ea. (l22|l for 
a fixed number of atoms A'^ = 10, 20, as function of the pump- 
ing 1/rj (in units of 1/k). The solid (dashed) line indicates 
the stable (unstable) solutions. The dimensionless parame- 
ters are Uo = 0.62 and Ac = 5, and the sites of the cavity 
potential are K — 50. The region shaded in grey corresponds 
to the values n < 7 where the tight-binding and single-band 
approximations are not valid. 



D. Validity of approximations and typical 
experimental parameters 

Various approximations have been imposed through- 
out the paper and it is in order to summarize them and 
discuss their reliabilities. The very first assumption of 
our model system is the adiabatic elimination of the ex- 
cited atomic levels and neglecting spontaneous emission. 
The justification of leaving out spontaneous emission, 
and at the same time taking into account for the cav- 
ity spectral line width and the single mode assumption, 
was discussed in detail in There we found effective 
atom-field couplings 50 ~ Stt x 100 MHz or go ^ 2tt x 700 
MHz using, respectively, the cavity decay rates k of [l3| 
or |18|]. These are presently slightly outside the regime 
of experimental reach. This may, however, be circum- 
vented by increasing the total number of atoms N [2^ . 
To assure the elimination of the excited atomic state we 
should have Aq ~ lOgov^) where the average number of 
photons n is typically 100. This suggests Uq < go and in 
our numerical analysis Uo is of the same order as k, while 
for realistic parameters one has go > k and thus ensure 
that the upper atomic state can be eliminated. 

We further imply the Gaussian, single band and tight- 
binding approximations in order to derive the many-body 
Hamiltonian ^ . The validity of these approximations is 
thoroughly analyzed in [l^l where it is found that they 
all break down for ?/ > 1, or equivalently \V\ < 1. In this 
work we consider y < O.b (or otherwise explicitly pointed 
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out) to meet this constrain. Especially, in the instabil- 
ity of field intensities as the pump amplitude is varied, 
shown in Fig.[3l are these approximations assumed to be 
violated. However, instability is still predicted, jumping 
from a strong cavity field to a weak one even if the true 
full state of the system cannot be predicted. Related is 
the conjecture of the Fermi sea state of the atoms. This 
requires an adiabatic change of any system parameters, 
but we have verified that our results are reproducible also 
for deviations from the Fermi sea state and therefore not 
restricted to the fully adiabatic limit. Consequently, the 
system is also robust against small but nonzero temper- 
atures. 

In solving for the dynamics, the set of Heisenberg equa- 
tions of motion (fT3|) and (fT4|) is truncated to contain the 
equations for a, and fi, which is believed to be justified 
in the regime of fairly large cavity decay considered here. 
The analytical steady state result for the field amplitude 
uses [a,g{n)] « dg{h)a and also neglects small terms in 
the final expression ([221) • This is supported since the 
results of the steady state numerical calculations of the 
Heisenberg equations coincide well with ones predicted 
by (|22p . even for small field amplitudes n < 10. 

To get an idea about characteristic experimental scales 
and quantities, we consider potassium atoms '^'^K and 
a wave length A = 800 nm, giving the recoil frequency 
Wr « 50 kHz. Thus, t = 1 corresponds to an unscalled 
time t = 20 fis, and estimating the critical temperature 
as T = hujr/kB one finds T ^ 0.4 fjK. Note that cavity 
decay rates used in this paper are of order of the recoil 
frequency, i.e. realistic but rather small. However, we 
have verified that k can be increased, while Uq, i] and Ac 
are scaled accordingly, and our results still remains. 

IV. DISCUSSION AND CONCLUSIONS 

We have reported on instability in the intracavity mean 
photon number, which are driven by the pump intensity 
and by the atomic density. Differing from usual opti- 
cal bistability, these nonlinear effects arise in the regime 
in which the atomic polarization is linear in the inten- 
sity of the cavity field. In fact, such instability is due to 
nonlinear coupling between the atomic motion and the 
cavity field, whose dynamics and steady state depend in 
a complex way on various parameters, such as pump am- 



plitude, frequency, and the atomic density distribution. 
The system we here consider is a gas of ultracold Fermi 
atoms, and the nonlinearity we observe is exclusively due 
to the coupling of the atoms with the resonator, whereby 
particle-particle collisions are neglected due to the Pauli 
principle. The atomic density distribution hence enters 
in determining the strength of the coupling of the atomic 
fermions to the cavity field, accordingly the height of the 
cavity potential which finally, closing this nonlinear cir- 
cle, determines the atomic density distribution itself. 

The dependence of the cavity field dynamics on the 
atomic density is indeed peculiar: At different atomic 
densities, different behaviours and instabilities of the cav- 
ity field are observed. One could hence conjecture that 
this system may also exhibit fiuctuations and instabili- 
ties in the mean number of atoms, in the regime in which 
this is not fixed, as it could be observed when the gas is 
put in contact with a particle reservoir. 

Further interesting outlooks emerge when considering 
the situation at larger intracavity intensities, in which the 
atomic polarization is nonlinear in the field. In this case, 
this nonlinearity, which corresponds to the one of disper- 
sive optical bistability, adds to the one of the coupling to 
the motion, and gives rise to a competition between cav- 
ity potentials. In fact, at lower intensities the potential is 
of the form cos^(A'), while at larger fields it is a superpo- 
sition of cos^(A') and cos^(X). Hence, instabilities in the 
atomic ground states may emerge, which are expected to 
give rise to novel behaviours. Such rich scenario is indeed 
exciting and requires different approaches to the study of 
this problem, which is under current investigation [29|. 
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